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> ■ Abstract 
O 

Q ! One-dimensional fractional statistics is studied using the Calogero-Sutherland 

\ model (CSM) which describes a system of non-relativistic quantum particles 

interacting with inverse-square two-body potential on a ring. The inverse- 
square exchange can be regarded as a pure statistical interaction and this 

^ , system can be mapped to an ideal gas obeying the fractional exclusion and 

O 

^ ' exchange statistics. The details of the exact calculations of the dynamical 

^ \ correlation functions for this ideal system is presented in this paper. An 

, effective low-energy one-dimensional "anyon" model is constructed; and its 

correlation functions are found to be in agreement with those in the CSM; 
and this agreement provides an evidence for the equivalence of the first- and 
the second-quantized construction of the ID anyon model at least in the long 
wave-length limit. Furthermore, the finite-size scaling applicable to the con- 
formally invariant systems is used to obtain the complete set of correlation 

exponents for the CSM. 
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I. INTRODUCTION 



The fractional statistics in low-dimensional (< 3) quantum systems had been a subject 
appreciated only by a few [0,0; however, the discovery of the quantum Hall effect has 
perhaps changed the status of the subject forever, promoting and reshaping it into one of the 
most profound as well as popular field in modern physics and into a necessary conceptual 
tool in condensed matter physics. 

The fractional statistics is usually discussed in the context of two-dimensional systems 
where the adiabatic transport of a test particle around the others can be used to deter- 
mine the statistics independent of the dynamical nature of the interacting system. In one- 
dimension (ID), however, the dynamics and the kinematics can not be decoupled in an 
unambiguous way (i.e., an exchange necessarily involves a scattering), and the assignment 
of the statistics to ID particles is largely a matter of taste. On the other hand, a good taste 
will yield fruitful concepts and tools. 

In this paper I study the ID fractional statistics using a specific model called Calogero- 
Sutherland Model (GSM) which describes a systems of non-relativistic quantum particles 
interacting with a pairwise potential that falls off as inverse- square of the distance between 
the particles. One of the most important practical features of this model is that the inverse- 
square potential can be regarded as a pure statistical interaction and the model maps to 
an ideal gas of particles obeying the fractional statistics Usually, the first step in 

understanding a general class of interacting Fermi system known as Fermi liquid is to study 
the ideal Fermi gas which gives rise to important concepts and tools such as the Fermi surface. 
Much in the same spirit I start with the simplest ID system (i.e., the GSM) obeying the 
fraction statistics. 

First, I list here some of the recent developments. A lattice cousin of the GSM called 
Haldane-Shastry Model (HSM), which corresponds to the SU(2) Heisenberg spin chain with 
the inverse-square instead of the usual nearest neighbor exchange, has triggered a surge of 
interest in this class of models 0. The HSM is known to possess the quantum symmetry 
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algebra known as Yangian |^ and can be considered as a model of ideal SU{2) spinon gas 
obeying the semionic fractional statistics 0. The SU{n) versions of the CSM |]T0[ and the 
HSM p!0| . [Tl[] now exist and in particular the spectrum of SU {n) HSM is known to possess 



the Bethe ansatz string structure |T2[. The list goes on, but I concentrate on the U{1) CSM 
in this paper. 

The CSM is intimately related to the circular ensembles in random matrix theory first 
introduced by Freeman Dyson |jl3|. In particular, the eigenvalue distribution functions for 
the orthogonal, unitary, and symplectic random matrices correspond to the ground state 
wavefunctions of the CSM at the interaction parameters A = 1/2, 1 and 2, respectively. 
(A = 1 case corresponds to the free Fermi gas.) Some static correlation functions of the 



CSM can be calculated using the techniques developed for the random matrices JM . 

More recently, Simons et. al. have been successful in mapping the CSM to the matrix 
model and to the non-linear sigma model where the supersymmetric algebra is applicable 
p!5| , and therefore are able to calculate the dynamical density-density correlation functions 
(DDDCF) for the CSM at A = 1/2, 1 and 2 [|T^. Haldane and Zirnbauer using the similar 
method calculate the one-particle Green's function at A = 2 (i.e., the symplectic case) 
P^ . Forrester also calculates some static correlation functions at integer interaction 
parameters using a generalized form of the celebrated Selberg integral formula [l^ . 



It turns out that the eigenstates of the CSM can be written in terms of Jack polynomials 
PJ20[| whose known algebraic properties provide a powerful and direct method for calculating 



the most general correlation functions. Recently, the author has been successful in calcu- 
lating the exact DDDCF and one-particle Green's function at arbitrary rational interaction 
parameters [^. The method employed is new and is one of the main subjects in this paper. 

The exact calculations of the correlation functions further provide conclusive evidences 
of the inherent fractional exclusion and exchange statistics embodied in the CSM 0. I 
also make a direct connection of this model to the edge states of the fractional quantum 
Hall droplet by first constructing an effective low-energy one- dimensional "anyon" fluid 
model based on a general gauge invariance argument, and then showing that the dynamical 



correlation functions agree with those of the CSM. This connection between the CSM and 
the edge states has previously been suggested by various people [ pi}] . 

This paper is organized in the following way. In Section y I show how to construct 
general eigenstates of the CSM, and in Section ^TT| some properties of the Jack polynomials 
are introduced. In Section |V| some key aspects of the exclusion and the exchange statistics 
are discussed. I show how to calculate the DDDCF and the one-particle Green's function 
in Section and |V|. In Section [VI I| the "harmonic-fluid" description of the anyon fluid is 



constructed and the dynamical correlation functions calculated; and I further show that they 
agree with those of the CSM. This agreement provides an explicit connection between the 
CSM and the system of coupled left- and right-edges of the fractional quantum Hall effect 
and further shows the equivalence between the first- and the second-quantized construction 



of one-dimensional "anyon" gas at least in the long- wavelength limit. In Section [VII1| the 
finite-size scaling as applied to the conformally invariant systems is used to obtain the 
complete set of correlation exponents. I also discuss some aspect of the lattice cousins of 
the CSM in Section HXI. 



II. EIGENSTATES OF CALOGERO-SUTHERLAND MODEL 



In this section I introduce the CSM and show how to construct the general eigenstates 



of the model following Sutherland |22|. First, the Hamiltonian for the CSM on a ring of 
length L is given by 

H_ J2 ^\ I l^^^'^K (2.1) 
where Tl^ /2m = 1 and d{x) = (L/27r)| sin(7rx/L)|. The ground state wavefunctions of the 



model at integer values of A corresponds to ID versions of Laughlin's wavefunctions p3[ and 
are given by 

^o = n(^«-^.)'n^fe"' (2-2) 

i<j k 
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where the current Jo = —X{N — l)/2 and zj = exp{i2'n'Xj/L). When < A < 1 there is 
another possible ground state with power 1 — A; however, only the solution with power A 
will be considered for reasons of continuity with A > 1 solutions. 

The excited states of this model are constructed by multiplying some symmetric poly- 
nomials to the ground state wavefunction, and this construction is analogous to that of the 
gapless edge excitations of the quantum Hall effect p4|. A general excited state \E'^ = ^^o^n 
is labeled by the quantum numbers n = (ni, n2, . . . , nAr), and satisfies the following new 
eigenvalue equation 

HJ^ = EnJn: (2.3) 

where H = Hq + XHi, and 

TV 

^o = E(^A)'' (2-4) 
i=i 

H, = Y.'f±^{zA,-zAJ. (2.5) 

j<k 

The eigenstates of the new Hamiltonian H are represented in terms of the following bosonic 
basis states 

N 

<^'(n) = En^r' (2-6) 

p j=l 

where the sum extends over all permutations of the integer set n which can be considered 
as a set of bosonic quantum numbers with no restrictions on their values. Since $(n) does 
not depend on the ordering of the quantum numbers, let ni > 77.2 >...> njv without loss 
of generality. These symmetric polynomials form a complete basis. 
The action of H on $(n) can be easily calculated and are given by 

HoHn)=^n'}j^n), (2.7) 

$(n) + 2 E ^(•••,n,-s,...,nfc + s,...) . (2.8) 

s=l J 

Hq generates only the basis state $(n) itself while Hi is responsible for generating a family 
of states $(. . . , — s, . . . , nfc + s, . . .) which are obtained from $(n) by all possible pairwise 



"squeezing" of the quantum numbers. If a state is generated from another by squeezing a 
pair of quantum numbers by one unit (i.e., rij —>■ rij — 1, rik —>■ rik + 1 for rij — Uk > 2), then 
I call the former "daughter-state" and the latter its "mother-state." 

The family of states can be organized into levels such that the members of a given level 
are mutually not related or unreachable (i.e., no mother-daughter relationship exists between 
the members in the same level) and the daughters of a member from a given level always 
belong to a lower level in the family. The highest-level mother-state denoted by |z/)i, where 
1/ is a level index and is equal to the total number of levels in the family, generates the 
entire family of daughter-states which are denoted by \fi)a where 1 < fi < u and a an index 
for the states in the fith level. While the order of the levels in the family can be uniquely 
determined, the order within a level is quite arbitrary. 

To provide an illustration of the above mentioned family structure, I give the following 
example. Let the highest-level mother-state be |6)i = $(4,3,1,0). Then, the members of 
the family are |5)i = $(4, 2, 2, 0), |4)i = $(4, 2, 1, 1), |4)2 = $(3, 3, 2, 0), |3)i = $(3, 3, 1, 1), 
|2)i = $(3, 2, 2, 1), and |l)i = $(2, 2, 2, 2). A pictorial representation of the family structure 
is shown in Fig. Q and I call it a level diagram. In Fig. |I| the states are represented by 
dots and each arrow connects a state A to B where B is reachable from A by squeezing on 
a pair of quantum numbers for A by one unit (i.e., connects a mother and her daughter). 
The levels are ordered from top to bottom, from the most squeezable to unsqueezable states, 
such that the highest-level mother-state is at the top; and therefore the arrows always point 
downwards and never upwards. A set of arrows that is topologically equivalent to a directed 
line forms a path. A weight W of an arrow is defined to be m{ni)m{nj) , where rij and Uj 
are the two quantum numbers squeezed to produce a daughter and m(/) the multiplicity of 
I in the quantum number set specifying the mother-state. In Fig. |I] the numbers next to 
the arrows are the corresponding weights. The states in the same level are not connected; 
hence, the Hamiltonian is diagonal in that subspace. 

A subfamily of a family can also be constructed by starting from a given state, which 
would be the highest-level mother-state of that subfamily, and grouping all her reachable 
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off-springs. The total number of subfamilies is therefore equal to the dimension (i.e., the 
total number of basis states) of the family. 

The lowest-level states (the unsqueezable states) clearly have to be of one of the two types: 
(I) $(. . . ,m,m,m,m,m, . . .) or (II) $(. ..,m,m,m,m,m — l,m — l,m — l,m — l,...). The 
type-I state corresponds to the ground state up to a global Galilean boost and the type- 
II to states with a single hole excitation. These one-hole states are eigenstates of H and, 
furthermore, since the following superposed state 

N 

^{x) = l[i^-z,), (2.9) 

i=i 

where z = exp(i27ix/L), can be expanded purely in terms of the one-hole eigenstates, \&(a;) 
describes a state with a hole localized at x. 

The matrix representation of H acting on the partially ordered state space is always 
triangular since the action of if on a given state always generates states belonging to lower 
levels. The eigenvalues, therefore, are simply given by the diagonal matrix elements. In 
particular the energy of an eigenstate spanned by a family with the highest-level mother- 
state <l>(n) is given by 

N 

j=l j<k 

The off-diagonal elements are given by 

, ,~, , I [Ep(n..P w.)]i^;« if. >/., 

a{li\H\v)p=l (2.11) 

I if < /i, 

where E^^ = '2XJ2j<k{^j ~ ^k), and the sum is over all possible path V from |z/)^ to \fi)a 
and the product over all the weights Wj of the intermediate arrows belonging to V. 
For the example given in Fig. |l|, H is represented by 
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H 



(1,1) ^(1,1) ^(1,1) ^(1,1) ^(1,1) ^ 

^(1,1) ^(2,1) ^(3,1) ^(4,1) ^(4,2) ^ 

n. A2^) (2,1) (2,1) (2,1) 

^ ^(2,1) ^(3,1) ^(4,1) ^(4,2) ^ 
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/ 



(2.12) 



n(i) 



i{Ijl\H\v) For example, e^g'-^j = 7 x E^^^^-^^-^ 



since there are five different 
ways to get from |6)i to |3)i with the following corresponding weights (see Fig. |l]): (a) 
|6)i - |4)i ^ |3)i; (b) |6)i ^ |5)i ^ |4)i ^ |3)i; (c) |6)i ^ |3)i; (d)|6)i ^ 14)^ ^ |3)i; (e) 
|6)i-|5)i^|4)2^|3)i. 

The eigenenergy given by Eq. ( |2.10| ) plus the ground state energy can be rewritten in 



terms of newly defined pseudomomenta kj as 

,.2 N 

E 



— VP 

2m ^tl ^ ' 



where 



Lkj = 2nlj + 7r(A — 1) X! ^S^i^j ~ ^i)- 



(2.13) 



(2.14) 



1=1 



The quantum numbers Ij are now distinct (half-odd) integers and are related to n/s by 
Ij = nj + {N + l-2j)/2. 

The distribution of kj determined by Eq. ( |2.14| ) is used to construct pictorial representa- 
tions of the eigenstates called motifs which are crucial for exposing the fractional statistics 
obeyed by the elementary excitations of the model. Detailed discussion of this subject is 
given in Section |IV|. 

III. JACK SYMMETRIC POLYNOMIALS 



The polynomial solutions of Eq. ( p^.3D in Section |I| is also known in mathematical 
literature as Jack polynomials ||25|. In fact, Stanley 12^] has shown that the complete set of 
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linearly independent solutions of Eq. ( |2.3| ) is indeed given by Jack polynomials up to global 
Galilean boosts (i.e., up to the factor Ojii ^/ where J is the current and takes an arbitrary 
real number) 0. 

For better readability this section is divided into two subsections: the first introduces the 
conventional notations used in mathematical literatures and the second some of the general 
properties of Jack polynomials. 



A. Introduction to notations 

Partitions are defined as sequences of non-negative integers in non- increasing order and 
are used to label the symmetric polynomials. They are denoted by bold-face Greek letters 
as 

K = (fi:i,K2, ...,KAr), (3.1) 

where «i > /t2 > • • • > i^n- Non-zero hij are called parts of n whose length (i.e., the total 
number of non- zero parts) is denoted by £(k). The weight of the partition is defined by 
1*^1 = Yll^=i i^j- If /ti + . . . + > /ii + . . . + /ij for alH > 1, then k > /x. 

Young diagram V^k) is used to graphically represent a partition: V^n) = : 1 < 

i < 1 < j < Hi}. The cell labeled by {i,j) is situated in the i-th row and the j-th 

column of the Young diagram. The diagram of k, therefore, consists of i{K,) rows of lengths 

Kj. 

A conjugate of n is denoted by ) and corresponds to a partition whose 

diagram is obtained by changing all the rows of T>{k) to columns in non- increasing order 
from the left to right. For example, the conjugate of k = (5,2,2, 1) is k' = (4,3, 1, 1, 1). 
Now, the following simple but useful identity can be derived 



In order prove Eq. ( |3.2| ) every cell in the ith row of I^(k) is filled in with an integer i — 1. 
Since n{K,) corresponds to the sum of all the integers in the diagram, the two different 
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expressions for n{K) are obtained depending on whether the numbers in each row or column 
are summed first. 

For a given cell s = of a diagram T>{k,) there are corresponding arm-length a{s) = 
Ki — j, arm-colength a'{s) = j — 1, leg-length l{s) = k'j — i, and leg-colength l'{s) = i — 1. 
The upper and lower hook-lengths are defined, respectively, as 



The symmetric polynomials are indexed by the partitions. The bosonic basis functions 
$(n) of the CSM are called the monomial symmetric functions and the quantum numbers n 
correspond to the partitions defined in the previous subsection. Since the quantum numbers 
are allowed to be negative integers, the correspondence is only up to some trivial constant 
translation or the global Galilean boost. This restriction to non- negative integer parts is 
more of a convenience than a restriction since the CSM Hamiltonian is invariant under the 
global Galilean transformation. 

I shall denote Jack symmetric polynomials as jy^(zi, . . . , z^) which are the solutions of 
Eq. (|2.3| ). If A = 1, Jack polynomials reduce to Shur functions which describe the excited 
states of the free fermions. At A = 0, it becomes the monomial symmetric function which 
is just the free bosonic wavefunction. If A = 2 or 1/2, they are called the zonal spherical 
functions. As A ^ oo, J^-^ reduce to the elementary symmetric functions. 

One way of defining Jack polynomials is through the differential equation ( |2.3| ). The 
other is based on the properties of the power-sum symmetric function p^, = Pk.iPk2'Pkz ■ ■ ■, 
where = J2j z'^" ■ Define a bilinear scalar product on the vector space of all symmetric 
functions of finite degree as 




(3.3) 



(3.4) 



B. General properties of Jack polynomials 



(3.5) 
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where zi^ = Y[i>i i^^^rriil, and = mi{n) is the number of parts of k equal to i. Using this 
definition, Macdonald |^8[ proved that there are unique symmetric functions satisfying the 
following three properties: 

1. Orthogonality: {Jk, J^)i/x = where is the normalization constant. 

2. Triangularity: Ji^ = J^fj, vt^^^^^i), where Vf^n = unless k < /j,. 

3. Normalization: If = n, then ^^,^1^ = n\, where k = (1, 1, ... , 1). 

n 

Jk are, then, constructed by Gram-Schmidt orthogonalization relative to the scalar product 
on the ring of polynomials. Stanley pH] proved that the normalization constant is given by 



Jk = n hUs)h^{s). (3.6) 
seK, 

There is another scalar product on which Jack polynomials are orthogonal: 



(k|/x)i/A = ( n / I jH^i^l^ ^2, • • • , ZN)j]i^{Zi, Z2, • • • , Zn) \{ \Zi - Zj 



2X 



where Zj = exp{i2Trxj/L) and A% = (1/L)^r^(l + A)/r(l + \N), and the bar over the 
polynomial denotes the complex conjugation. Note also that the multidimensional integral 
above is equal to times the constant term in 

J'J\l/z,, l/z„ . . . , l/z^)4/\z„ z,,..., z^) n(l - -)^ (3.8) 



Eq. (p.7|) has been conjectured first by Macdonald and then later proved by himself ||28 
and also by Kadell |^ . 



Since Jack polynomials span the vector space of symmetric functions, they can be used 
to expand any symmetric functions. This property is particularly useful in calculating the 
correlation functions of the GSM as will be shown later in this paper. Here are some of them 
PT|,B1: 
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E^r = y E ^^JU\z„...,z^), (3.9) 
i=l ^ \n\=n Ji^ 

11(1 - z,r = E W^^^/^^i, • • • , zn), (3.10) 
j=i n ^ 

where [a]^ = n(^,,)g^v{a + (j - - (« - 1)}, and {a}^ = n(i,,)eAi{a - - 1) + (j - !)}• 
The sum in Eq. ( ^.1(J| ) extends over all possible partitions while in Eq. ( |3.y| ) it is restricted 
to partitions with weight \k\ = n. The prime in [0']^ denotes that the product does not 
include the cell (0, 0); otherwise the total product is trivially equal to zero. also satisfies 

JII^{wZi,WZ2, . . .,wzn) = w^'^^jII^{zi,Z2, . . .,zn), (3.11) 
since Jack polynomials are homogeneous functions of degree |k;|. 



IV. FRACTIONAL STATISTICS 

I divide this section into two subsections. In the first (second) subsection the fractional 
exchange {exclusion) statistics is discussed in the context of one-dimensional models. The 
exchange statistics in two-dimension is directly relevant to the fractional quantum Hall effect 
and various people have made contributions to this fascinating field p3. In one- dimension, 



however, the definition of fractional exchange statistics is rather obscure and incomplete with 
a possible exception of the CSM. One the other hand, the definition of fractional exclusion 
statistics is spatial dimension independent and is based on the structure of the Hilbert space 
rather than the configuration space. While the fermions obey the well-known Pauli exclusion 
principle, more exotic particles may obey a "generalized exclusion principle." p[ 



A. Exchange Statistics 

The first full mathematical treatment of the fractional statistics is due mainly to Leinaas 
and Myrheim who used the multiply connected topological structure of the configuration 
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space of collections of identical particles to show the possibility of exotic statistics in spa- 
tial dimension less than three 0. While their idea has been extensively applied to two- 
dimensional systems especially in the context of the fractional quantum Hall effect, very 
little attention has been paid to the one-dimensional (ID) systems. Perhaps the main dif- 
ficulty in ID systems is that the physical exchanges of the particles necessarily involve 
scattering processes and that there is no known unique way of un-tangling the kinematic 
aspect of the fractional statistics from the dynamical processes. 

For integrable ID quantum systems, however, there is a general consistency condition 
known as the Yang-Baxter equation (YBE) which essentially puts strong restrictions on the 
scattering matrices. Therefore, intuitively the braiding of the particle "trajectories" in one- 
dimension when properly defined may be given by the YBE . In fact the YBE is known 



to be intimately related to the the knot theory and the braid groups |^5[. Hence, perhaps 
the following quantum Yang-Baxter equation should be interpreted as the one-dimensional 
generalized braiding relations. 



Suiv - u)Si3{v)S23{u) = S23{u)Si3{v)Si2{v - u) , (4.1) 



where Sij{v) is the scattering matrix in the tensor product of linear vector spaces, V^V^V 
and acts non-trivially only in the ith and jth space, e.g. Si2{v) = R{v) I where R{v) 
is a matrix defined in V ® V. The parameters v and u are called spectral parameters and 
are equal to the rapidity differences between two colliding particles. One can rewrite Eq. 
( [4. ID in terms of R{v) = PR{v) where P denotes the transposition, Px ®y = y®x. Define 
matrices Ri{v) = I®---®R{v)®---®IonV®---®V, where R{v) acts on i-th and i + 1-th 
spaces. The matrices Riiy) satisfy the following relations 

Ri{v)Rj{u) = Rj{u)Ri{v) if|i-j|>2, (4.2) 
R^+i{v - u)Ri{v)Ri+i{u) = Ri{u)Ri+i{v)Ri{v - u). (4.3) 

Without the spectral parameters the R matrices satisfy the braiding relations. The complete 
description of the fractional statistics based on the YBE needs to be worked out. 
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At present it seems unclear to me how the construction based on the YBE can be apphed 
to a general ID (i.e., not necessarily integrable) systems. Maybe there is a way to map a class 
of general ID models to integrable models with in some sense small residual interactions that 
break the integrability. Instead of pursuing this incomplete Yang-Baxter story any further, I 
give here somewhat heuristic but intuitive device for constructing the ID fractional statistics. 
The A^-particle state is constructed as follow 

dxi ■ ■ ■ dxN'^ixi, . . . , xn, t)(p*{xN, t) ■ ■ ■ 0*(xi, t) |0), (4.4) 

where (j)*{x,t) is the anyon creation operator. The "anyon" fields satisfy the following equal 
time commutation relations 

[0(x, t), (Piy, t)]x = [<P*{x, t), <l)*{y, t)]x = [0(x, t), <p*{y, t)], = Oif x ^ y, (4.5) 

where [x, y]x = xy + exp(?7rAsgn(?/ — x))yx. 

The wavefunction \I'(xi, . . . , xat, t) is multi-valued and satisfies 

^(. . . , x,, . . . , X, . . . , t) = e-^sgn(x,-x,)^^_ . . , x„ . . . , x,, . . . , t). (4.6) 
The integrand in Eq. ([4. 41), however, is always single valued since the phases arising from 



the wavefunction and the "anyon" field operators are set to cancel each other. 

If I specialize to the CSM the fractional exchange statistics can be formulated in the first 



quantized language. First, the Hamiltonian (|2.1| ) can be rewritten as 

2 



1 ^ / .vr/iA ^ 



niXj - Xk 



P,k\ , (4.7) 



where pj = —ihdx^ is the momentum operator and Pij the particle exchange operator. The 
extra term added to the momentum operator is an ID analog of the Chern-Simons gauge 
field. 

In two-dimension there are two well-known ways to code the fractional statistics for the 
ideal anyon gas. One is to take the free Hamiltonian and require that its wavefunctions be 
multi- valued as in Eq. ([4.6|) . The other is to introduce the Chern-Simons gauge field and 
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write the Hamiltonian in terms of this gauge field and further require that the wavefunctions 
be single valued and symmetric. The Hamiltonian ( ^.Tp corresponds to the ID version of 
the second formulation of the fractional exchange statistics. Hence, the following symmetric 
wavefunctions are the eigenstates of the ID anyon system 

^U^i, ...,xn) = <f{xi, xm)'^oJL (4.8) 

where the "ordering function" (p is introduced to maintain the total wavefunction symmetric 
In particular ip keeps track of the braiding of the particles and is set to cancel the 
exchange phases arising from \l/o. I require the ordering function to satisfy 

p{xi, . . . ,Xm,- ■ ■ ,xn) = e*('""^)''^v9(xi, . . . , . . . , Xtv), (4.9) 

where the hat over the variable denotes the absence of that variable from the function. In 
other words, to remove the particle at Xm it is necessary to pass through m — 1 particles. 

In Section |VII| I construct explicit second quantized "anyon" fields and show their con- 
sistency with the first quantized formulation of the fractional exchange statistics for the 
CSM. 

B. Exclusion Statistics 

The notion of fractional exclusion statistics based on the so called "generalized Pauli 
exclusion principle" has first been formulated by Haldane and applied to the elementary 
topological excitations of general condensed matter systems 0. This new concept of statis- 
tics is based on the structure of the single particle Hilbert space of the elementary excitations. 
More specifically, the change in the size of the available states (AD) in the Hilbert space 
as the number of particles (i.e., the elementary excitations) is changed (AA^) for a given 
system with fixed boundary condition defines the statistics of the particles with the statisti- 
cal parameter defined as g = —AD/ AN. Hence, for example, the bosons and fermions are 
identified with g = and g = I, respectively. 



15 



In order to facilitate proper understanding of the exclusion statistics in the context of the 
CSM I introduce a pictorial representation of the eigenstates and make the identification of 
the excitation contents of the states easier 0. Eq. ( p.l4[ ) gives the occupation configurations 



of the pseudomomenta kj for all the eigenstates of the CSM. The quantum numbers {Ij} in 
Eq. ( [^.141 ) are distinct (half-odd) integers and in the ground state are given by the following 

set 

{.,.(-^,-^,...,^,^1 ,4.10) 



Therefore, the ground state pseudomomenta are given by 



{*«} = -^(A^ - 1). -^(iV - 3). , ... ^(A' - 3), ^(A^ - 1) . (4.11) 



The total ground state energy E° = Eji^jY is equal to tt'^X'^N{N^ - 1)/3L^. The excited 
states are given by integer displacements of {Ij}- Therefore, two neighboring pseudomo- 
menta for any arbitrary state must be separated by 

Akj = \k^- kj_i\ = ^(A + 0, (4.12) 

where / is a non-negative integer. 

In order to construct a picture that exposes the excitation content of the excited states, 
I let A to be a rational number p/q with p and q coprimes and introduce one-dimensional 
lattice with the lattice spacing equal to 27r/gL. I assign each lattice point with 1 if that 
lattice point coincides with the value of one of the occupied pseudomomenta and with 
if it does not. Hence, the ground state for A = 3/2 and = 10 is represented by 
. . . 0000000010010010010010010010010010010000000 . . .. All the other excited states can be 
obtained from this ground state configuration by displacing the ones such that the number 
of zeroes between any pair of ones is equal to p — 1 + ql where / is a non-negative integer. 

I use three different names for the particles in the model — real, pseudo, and quasipar- 
ticles. The real particles are, of course, the physical quantum particles described by the 
canonically conjugate coordinate and momentum variables {xj,pj}. The pseudoparticles 
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are described by the pseudomomentum operator (see Eq. ([4.7|)) which a sum of the usual 
momentum and the ID "statistical gauge field." The quasiparticles are the elementary exci- 
tations of the system. Because the pseudoparticles form an ideal gas, the quasiparticles are 
essentially same as the pseudoparticles excited out of the condensate. The holes left behind 
in the pseudoparticle condensate will be called quasiholes. Hence, the name "pseudo" and 
"quasi" will be used interchangeably in some cases p. 

Since p — 1 zeroes are always required between the ones I call them bound zeroes which 
seem to represent the mutual statistical exclusions. In the case of the free Fermions carrying 
the fiux 2tt in units where e = h = c= lso that the fiux quanta $o = 27r, the minimum 
separation of the momenta is 2tt/L which can be considered as the mutual Pauli exclusion. 
In this case the minimum separation is 27rA; therefore, it is natural to assign 2ttX fiux 
attached to the pseudoparticles. The remaining zeroes are called unbound zeroes. The q 
consecutive unbound zeroes in the condensate of the pseudoparticles constitute a single hole 
excitation. Thus, if a pseudoparticle is removed from the ground state condensate, then 
there are p unbound zeroes in place where the one is removed. This state is forbidden if 
g 7^ 1. In general a minimum of q ones must be removed so that they leave behind at least 
pq unbound zeroes which break up into p holes. From the view point of the particles (holes) 
the change in the number of available single particle states is p {—q) while the change in 
the number of quasiparticles (quasiholes) in the system is —q (p). Therefore, the statistical 
parameter g for the quasiparticle (quasihole) is g = p/q = \ {g = q/p = 1/A). This inverse 
relationship between the statistical parameters of the particles and holes is essentially the 
Chern-Simons duality. To summarize the fractional exclusion statistics, q particle excitations 
are accompanied by p hole excitations. 

The configurations constructed above are representations of the diagrams of partitions 
©(k) introduced in Sec. PTT| . The part Hj corresponds to the displacement of jth quantum 
number from the ground state (i.e., Ij — Ij if Ii > I2 > . . . > In)- The excitations given by 
K, include only the states with non-negative displacements (i.e., kj moved only to the right) 
and all the other states are obtained by global Galilean transformations. Therefore, each 
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row (column) in the diagram corresponds to the particle (hole) excitations of the CSM. 
Following Yang and Yang [^] and Sutherland I can also construct the thermody- 



namics P,|38|. In the thermodynamic limit it is convenient to define the hole distribution 
Ph{k) and the particle distribution function pp{k). The function pp{k) are given by the solu- 
tions of Eq. ( [^.14D while Ph{k) by the corresponding complementary equations given by the 
unused quantum numbers /-,•; thus, they satisfy 1 = phik) + \pp{k). This equation states 
that one hole and 1/A particles (or A holes and one particle) have equal weight in occupying 
the volume in fc-space. Hence, for A 7^ 1 the particle-hole symmetry is broken; and in the 
bosonic case (A = 0) the symmetry is maximally broken. With a proper normalization one 
can equally state that 1 = {\/\)ph{k) + Pp{k). This particle-hole duality is the essence of 
the Chern-Simons duality. 

The thermodynamic function is given hj Q = E — TS — pN, where the energy (E), the 
particle number (N), and the entropy {S) are given by 

E = V J dkpp{k)e{k), (4.13) 
N = V J dkpp{k), (4.14) 
S = V j dk{{ph + pp)\og{ph + pp) - ph^ogph- pplogpp}. (4.15) 

Here, V is the volume of the system. By minimizing Q{pp,ph) with respect to the density 
functions the following relation can easily be obtained 

(1 - \pp{k))\l + (1 - X)pp{k)y~' = pp{k)e^<^^-^^l^. (4.16) 



This equation is also obtained by Wu using a different method [39 . 



V. DYNAMICAL DENSITY-DENSITY CORRELATION FUNCTION 

In this section I show that the dynamical density-density correlation function (DDDCF) 
can be calculated exactly for the CSM using the known properties of Jack polynomials. The 
DDDCF is defined by 
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;v(0|p(x',t')p(a:,t)|0)^ = ;v(0|e^^-*'p(xOe-^^-*'e^^-V(a:)e-'^-*|0)^ 

= ;v(0|p(x')e-^(^--^-)(*'-*V(a;)|0),v, (5.1) 

where the reducec? density operator p{x) = (l/L) J2j 6{x — Xj) — N/L and |0)Ar the normahzed 
A^-particle ground state. 

The first step in calculating the DDDCF is to expand p{x)\0)j\f in terms of the eigen- 
states of the CSM (i.e., Jack polynomials). The delta function 6{x) is a periodic function 
with a period L and thus can be expressed as a Fourier sum (l/L) X]m=-oo 6xp(227rxm/L). 
Therefore, I can write p{x) as follow 

-| oo 

P(^) = 7E(^>— + (5-2) 

^ m=l 

where z = exp(i27ra;/L), zj = exp{i2iTXj/ L), and Pm = Y.f=i z™"- The power sum can be 
expanded in terms of Jack polynomials using the identity Eq. (|3.9|) . 

Using the orthogonality relation Eq. ( p.7|) and its extension Eq. (|3.8|) , I obtain the 
following expression for the DDDCF 

;v(0|p(x, t)p(0, 0)|0);v = y^T^ E -^JW^TlV^ cos(2vr|/.|x/L)e-*^'^, (5.3) 

where E^, = (27r/L)2 ^^^^(k^ + A(A^+ 1 - 2j>j). The coefficient [0']^ in Eq. (U) vanishes 
unless the diagram T>{k,) has no more than p columns of length longer than q and q rows of 
length longer than p. In other words, the intermediate states contributing to the DDDCF 
has precisely p hole and q particle excitations. This is a conclusive evidence of the ideal 
fractional exclusion statistics the CSM quasiparticles and quasiholes obey. 

The DDDCF in the thermodynamic limit greatly simplifies as is shown in Appendix A; 
and it is given by 

(0|p(x,t)p(0,0)|0) = C f[ [j^dx^ n [j^dy^ F{q,p,\\{x,,y,}) cos{Qx)e-''^\ (5.4) 
where Q and the total momentum and energy, are given in units of fi and 'h? /2m by 

Q = 27rpo (e^. +E% | ' (5-5) 

Vi=i i=i / 
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E = (27rpo)' E + E ^H{y^j) . (5-6) 

Vi=i i=i / 

with po = ep{x) = x{x + A) and eniv) = A?/(l — y). Xj{ep) and Vj^eu) are normalized 

momentum (energy) of the quasiparticles and the quasiholes, respectively. The normaliza- 
tion constant C is given by 

r"(A)r"(i/A) ^/^ r(g-(j-i)/A) V 

^l^,n,Aj nr=ir^(P-A(.-i))n-=ir2(g-(,-i)/A)^iUr(i-(j-i)/A)J ' ^'-'^ 

y2p{q-l)Y2l \ 



Finally, the form factor F{q,p, X\{xi, yj}) is given by 



l/A 



The form factor has been conjectured by Haldane based on the clues given by the works 
of Simons et. al, Galilean invariance, and U{1) conformal field theory [Q. A less general 
form of the DDDCF at integer values of A has also been reported ||41|| . 

The region of support in the energy- momentum space for the DDDCF at A = 5/3 
corresponds to the shaded area in Fig. |. It is obtained by convoluting the dispersion 
relations of the q quasiparticles and p quasiholes as given by Eqs. (|5.5|) and ( |5.6| ) . At low- 
energies 2p distinct sectors indicated by darker shade emerge as expected of ID metallic 
system. (A generic ID system, however, will have an infinite number of these low-energy 
sectors.) 

It is also worth noting that there is a qualitative difference between the finite and the 
infinite system. The intermediate states represented by the diagrams with \K,j — Kk \ and {i^'j — 
n'j^l of order 0{N) are the only states that contribute to the DDDCF in the thermodynamic 
limit (i.e., iV, L — > oo with N/L fixed). I call this a super- select ion rule p. (There can 
be violations of this rule at the low-energy limit as discussed in Appendix A.) Perhaps, the 



phase transition that occurs in 2D QCD as — > oo is related to the CSM super-selection 
rule. 
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VI. ONE-PARTICLE GREEN'S FUNCTION 



Unlike the DDDCF calculated in the previous section, evaluation of the one particle 
Green's function depends on the actual statistics of the real particles since it involves the 
reordering of the particles. For example, the one-particle density matrix (the static limit of 
the Green's function) for the A = 1 case is trivially simple if the ground state \E'o is taken to 
be the fermionic wave function. On the other hand, if the wavefunction is taken to describe 
a bosonic system, the calculation gets quite complicated and different from its fermionic 



counterpart as shown in P3|. 



As described in section ^ I take the modified GSM given by Eq. (|4.7| ) as one-dimensional 



anyon system with the exchange phase given by exp(iA7r) and Eq. (|4.8| ) as the corresponding 
wavefunctions. If the mth particle is removed from the ground state of + 1 particles, the 
remaining term after factoring out the ground state wavefunction of the N particle system 
is 

( n - (Ui^m- zA ( n Z^'A Z-J^I\^^--^^-\ (6.1) 

\j{<m) j \i(>m) / \k+m j 

where the last term is determined from the property of the ordering function given in Eq. 



( [4.91 ). Therefore, up to an overall phase factor that does not depend on the position of the 
removed particle, the destruction operation \E'(x) on the ground state of + 1 particles is 
given by 

A ^ 

*(x)|0)^+i = - V^>)^, (6.2) 

where z = exp(z27rx/L) and A% = (1/L)^r^(l + A)/r(l + XN). The statistical phase 
arising from the ordering function makes the destruction operation symmetric with respect 
to the permutations of {zj}. The symmetric function njLi(2; — Zj)'^ can then be expanded 
in term of Jack polynomials using Eqs. ( |3.10| ) and ( p.ll|) . 



The hole propagator part of the one-particle Green's function is defined as 
^+l(O|vi/t(x^^0^(a^,^)|O)7v+l = iv+l(O|e^^-+^*'vi>t(x0e-^^-*'e^^-*vl/( 

= ^+i(0|^^(x')e-'«^^-^^)-^)(*'-*)^(x)|0)iv+i, (6.3) 
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where the chemical potential fi = -E'at+i — E^. By expanding \E'(a;) in terms of Jack polyno- 
mials and using the orthogonality relation (|3.7|) the propagator at finite N and L is evaluated 
to be 

A.+i(0|vI/t(x,t)vl/(0,0)|0)^.+i = (iV+ 1) (^)' 

where the additional + factor comes from the freedom of choosing one of available 
particles to destroy and create. The coefficient {— A}^ vanishes unless the diagram T'(k) 
has at most q — 1 rows of length greater than p and p columns of length greater than q — 1. 
Therefore, the intermediate states for the propagator is spanned by g — 1 quasiparticles and 
p quasiholes. This is a very important result. The exclusion statistics of the quasiparticles 
and quasiholes is completely consistent with the anyon statistics of the real particles. P| 

Taking the thermodynamic limit of Eq. (|6.4|) is almost identical to that of Eq. (|5.3| ) and 
is given by |^ 

m\x,t)^iO,0)\0) = PoDl[( rdx,)f[( [\) F(g-l,p,A|{x„y,})e^«Q-««)--(^-'^W, 

(6.5) 

where the chemical potential /i = (vrApo)^ and the back fiow Qq = vrApo- E{q—l,p, X\{xi, i/j}) 
is still given by Eq. ( |5.9| ) and D by 

A2p('?-i)r2(n) 

Q and E are same as before except for the number of Xj^s. At integer values of A (i.e., q = 1 



case where only quasiholes are excited), based on the equal-time results of Forrester fTS 
Haldane made a conjecture which agrees with this formula. The regions of support for 
the hole propagator is given by the shaded area in Q > (or Q < 0) in Fig. There are 
also shifts in E and Q by — /i and —Qo, respectively. 
It is also interesting to consider the following function 
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^i{x) = Y[{z-z,r, (6.7) 

i<j 

where m is a positive integer. This function can be expanded in terms of Jack polynomials 
using Eq. ( |3.10|) with coefficients containing the term {— wi}^ which vanishes unless the 



diagram T'(k;) has no more than m columns and no rows longer than m. Therefore, "^^{x) 
acting on the ground state creates exactly m quasiholes. "^^{x) is a generalization of the 
one-hole state given in Eq. (pl9|). The propagator can easily be calculated and its form factor 
is given by F{0, m, X\{yj}). This result is consistent with the conjecture that the minimal 
form factor for any two point correlation functions whose intermediate states contain only 
n quasiparticle and m quasiholes is given by F{n,m, X\{xi,yj}). 

VII. LOW-ENERGY EFFECTIVE THEORY AND COUPLED FRACTIONAL 

QUANTUM HALL EDGE STATES 



Following the standard method in Luttinger liquid theory |^ I construct an effective 
low-energy model for the ID anyon system. When A is an integer this system is equivalent 
to a coupled system of left- and right-moving edge states of the FQHE. The excitations 
on a single edge of the FQH fluid moves only in one direction because of the externally 
applied magnetic field; and they have been rather thoroughly studied using the so called 
chiral Luttinger liquid theory which is intrinsically anomalous When there is an 
extra edge on the FQH droplet (e.g. a strip, annulus, or cylindrical geometry instead of 
the disk geometry) and when the edges are close enough, new phase space opens up as a 
result of the fractional charge transfer between the edges. There actually have been many 



suggestions that the CSM is related to the edge states . In this section I show in another 
way that the CSM is an exactly solvable ID anyon system by calculating the correlation 
functions of the effective model and finding exact agreement with those of the CSM in the 
long- wavelength limit. 



A general gauge invariance argument |^ can be used to map out the qualitative structure 



of the excitation spectra in the energy-momentum space. I use the cylindrical geometry |^ 
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with coordinates x and y defined as shown in Fig. ^. I put some incompressible fluid perhaps 
made up of A anyons on the surface of the cyhnder with a confining potential in y-direction 
such that two edges are created along x-direction. (Of course, when A is not an integer the 
edges are not related to the FQHE edges which are more complicated composite type [0].) 
The cylinder with circumference L is also pierced by a thin solenoid along its longitudinal 
axis, thus inducing a magnetic fiux $ = AL through its cross-section and a vector potential 
A = Ax on its surface. There is also magnetic fields B normal to the entire surface of the 
cylinder. Since the left- and right-moving sectors of the GSM completely decouple in the 
low-energy limit, this geometrical construction is in fact equivalent to the GSM and is useful 
to show intuitively how the low-energy sectors are related. 

The fictitious fiux applied to the cylinder is used as a passive device for mapping out 
the regions of low-energy excitation sectors of the coupled anyon edges. By virtue of the 
gauge invariance the energy is degenerate at fiux equal to 27m (in units where h = e = 1) 
which corresponds to the momentum of 2TTnpo. The adiabatical change in fiux from to 27r 
induces an elementary excitation carrying charge —q/p and fiux —2tt (i.e., a quasihole) to 
move from one edge to the other; and, in between, it costs finite energy to the system since 
the bulk is incompressible. In fact the real space configuration of an incompressible droplet 



is equivalent to a "Fermi sea." [^] The edges of the fiuid coincide with the locus of points 
in which Fermi energy crosses the external confining potential. The only degrees of freedom 
left in "Fermi sea" in the low-energy limit is the "Fermi surface" fiuctuation. This "Fermi 
sea" obeys the fractional exclusion statistics. Since there is a symmetry along x-direction, 
the fiuid can be regarded as an ID "Fermi sea" and the left- and the right-edges as two 
"Fermi points." 

In light of the observations made so far it is reasonable to model the coupled edges by 
a full (non-chiral) ID anyon system. First, I construct the one- dimensional anyon creation 
operator as follow 

(x) = ^ya;)e'^^("\ (7.1) 
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where is the boson creation operator and A the statistical parameter. The operator 
giA6»(x) jg called "disorder operator" |Q that creates a kink (or vortex) of size ttA at 

position X and is defined as 



e{x) 



TT 



p{x')dx', 



(7.2) 



where p{x) is the density operator. Therefore, "^xi^x) is a composite operator that creates a 
boson plus a vortex; so, it is an anyon creation operator. 

The boson creation operator \E'^(x) for a system with multiple low-energy sectors has 



previously been constructed and is given by 



(7.3) 



m=—oo 



where the phase field 0(x) is defined by the following canonical commutation relation: 
[</)(x), p(x')] = i6{x — x'). The multi-sector density operator is also given by 



p{x) = p{x) E e 



i2me{x) 



(7.4) 



Here, the sectors are connected by the operator exp{i2m6{x)) which creates a vortex of size 
27rm. This dynamical device now replaces the passive device I previously used to map the 
system from one sector to the next by supplying external magnetic fiux to the cylinder. 

The fields 6{x), (f){x) and p{x) can now be expressed in terms of Tomonaga boson oper- 
ators b'^ and b as [^ , ^ 

1/2 



9{x) = 00 + Tcpox + ie'^ E 



7rk 



2L 



"'"■(&L, + bk) 



TC 



1/2 



p{x) = Po + e^J2 



k 



2kL 

1/2 



-''""ibl, - b,) 



27iL 



-ikx fji 



\bl, + b,). 



(7.5) 
(7.6) 

(7.7) 



I assume here that the most important interaction is the statistical interaction; so, I can set 
the Bogoliubov parameter e^*^ = 1/A. 

Now, (x) can easily be shown to satisfy the following anyon commutation relation 
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e'"^^g^("'-")^l(x')^i(x) for X ^ x'. 



(7.8) 



Eq. ( |7.8|) can thus be used as a defining relation for the one-dimensional anyons. 

The Hamiltonian is diagonal in Tomonaga boson operators and is given as usual by 



where Vg is the sound velocity. 

Using the operator identities (i) e^e^ = exp(e" — l)e^e^ if -B] = aB, (ii) e^e^ = 

exp{-[A,B])e^e^ and e^+^ = exp{-[A, B]/2)e^e^ if [[A,B],A] = [[A,B],B] = 0, the 
following correlation functions are calculated 



where = x^v^t, fi the chemical potential and the coefficients A^ (-Bm) are regulariz at ion- 
dependent constants. 

The Green's function (\l/^(x, t)\l/A(0, 0)) for the sector m = 0, where the charge transfer 
from one edge to the other is forbidden, is given by only the right-movers (or only the 
left-movers if were properly redefined) even though the anyon creation operator ^I^j(,(a;) 
contains both the right- and the left-moving bosonic modes (i.e., and bk for both A; > 
and A; < 0). This chiral sector emerges naturally in this theory without explicitly imposing 
the chirality condition. 

In an isolated chiral theory coupled to a gauge field the charge is in general not conserved 
and because of this the theory is known to be anomalous and not physical. As far as the 
isolated chiral theory is concerned the only physical sector is the chiral sector where no 
charge transfer between the edges is possible. One can, however, consider an "almost" 
chiral theory which describes a system of two chiral edge states that are independent except 



H = hvsJ2 \9\blbg, 



(7.9) 



q^O 




(7.10) 



(7.11) 
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for the charge leakage from one side to the other and vice versa. For the sake of definiteness 
I concentrate on the right edge and propose the "almost" chiral system with the following 
Hamiltonian and the field operators, 



H'' = nv,J2qblb, 



g>0 



e\x) = 00 + np^x - ^ E T 

1/2 



nk 



2L 



1/2 



k>0 



TT 



2kL 



^i2(m+A/2)0«(x) 



(7.12) 

(7.13) 

(7.14) 
(7.15) 



The chiral system is constructed only with the right-moving Tomonaga bosons. The Green's 
function in this case is given by 

^ ^ (m+A)/A 



(vI/f(a;,t)M/f(0,0)^p^ E Cr, 



J{2TTp^{m+\/2)x+ti^t) 



(7.16) 



The m = sector here is equivalent to the corresponding chiral sector of the non-chiral 
model. As expected only the right-movers contribute to the Green's function for all the 
other sectors. 

Asymptotic expansions of the correlation functions of the GSM have been calculated in 
Appendix B; and they agree with Eq. ( |7.1CI| ) and Eq. ( |7. 1 1| ) . This agreement between the 
correlation functions of the explicitly constructed anyon model and the GSM shows that 
the first- and the second-quantized construction of the anyons are in fact equivalent in the 
long- wavelength limit. 



VIII. FINITE-SIZE SCALING AND CORRELATION EXPONENTS 

There is another extremely elegant and powerful way of obtaining the exponents of 
the correlation functions in the long- wavelength limit. If the dispersion relations of the 
elementary excitations of the one-dimensional quantum model in the low-energy limit is 
linear (i.e., has Lorentz as well as Galilean invariance so that the space and time variables 
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are on equal footing), the principle of conformal invariance is applicable and, in particular, 
the finite-size corrections to the energy and momentum become universal and are directly 
related to the exponents of the correlation functions |^ . 



The finite-size scaling has previously been applied to the CSM by Kawakami and Yang 
1 5^ without the benefit of recently uncovered knowledge namely that due to the ideal 



fractional statistics the CSM particles obey the intermediate states for the DDDCF and the 
Green's function are spanned by finite number of the elementary excitations. The shaded 
regions in Fig. ^ are the relevant sectors for the correlation functions (for the Green's 
function only the positive or negative momentum sectors contribute). A complete set of 
relevant exponents is obtained using the finite-size scaling in this section. 

Let {pl, h^lpji, hfi) be a label for a low-energy sector spanned by pi {hi) left-moving and 
Pr {hn) right-moving quasiparticles (quasiholes). The relevant sectors for the DDDCF are 
(0, n|g,p — n) and {q,p — n\0, n) and for the Green's function (0, n\q,p — n) or {q,p — n|0, n) 
where n = 0,1,..., p. First, for the DDDCF the lowest-energy state in the sector m is 
characterized by m quasiholes on right (or left) side of the "Fermi sea" transfered to the left 
(or right). Hence, the pseudomomenta kj for this sector is given in terms of the ground-state 
pseudomomenta kj as kj = kj + 27rm/L; and the energy and momentum by 

E^ = E, + ^ — , (8.1) 
P„ = 2nmpo, (8.2) 

where the sound velocity Vg = 2ttXpq. The scaling dimension x = + h~ of the density 
operator and its conformal spin s = — h~ are therefore given by x = m?/\ and s = 0, 
where are right (left) conformal weights and 2h^ are the exponents that actually appear in 
the correlation functions for the right(left) movers. Thus, 2h^ = 2h~ = im? /\ in agreement 
with the exponents found in the previous section and in Appendix B. 

Now, for the hole propagator the energy and momentum of the sector m characterized 
by m quasihole transfers from the left to right "Fermi points" and thus the pseudomomenta 
kj = kj + 211771/ L, where j = 1, . . . , — 1, are 
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Em = Eo-fi + ^i — -n + -\ , (8.3) 
= 27r(n - A/2)po + ^ ~ ' ^^"^^ 

The chemical potential /i = (vrApo)^ is associated with the particle destruction. I remove the 
rightmost pseudoparticle with k% = 7tX{N — 1)/L from the ground state condensate and do 
not introduce any separate selection rules in contrast to Ref. |^2|. In this case the right and 



left conformal weights are different (due to the non-zero conformal spin s = A/2 — and 
are given by 2h^ = {m — A)^/A and 2h^ = m^/A as expected. The right and left conformal 
weights would be switched if the excitations were caused by the removal of the leftmost 
pseudoparticle at — 7rA(A^ — l)/L. 

Apparently, the long-range interaction does not destroy the conformal invariance of the 
CSM in the long- wavelength limit. This is expected from the linear dispersion relations of 
the CSM quasiparticles and quasiholes in this limit. 

IX. LATTICE MODELS 

The Haldane-Shastry model (HSM) corresponds to a lattice generalization of the CSM 
at A = 2 and has the following Hamiltonian 

H = JoEj}^y (9-1) 

where d{m) = (A^/7r)| sin('7rm/A^)| and Sj is the SU{2) spin operator acting on site j. The 
model possesses a quantum symmetry called Yangian and exhibits supermultiplets struc- 
tures whose spin contents are exactly reproducible from an asymptotic limit of the ther- 
modynamic Bethe-ansatz equations [1^. Furthermore, there is one-to-one correspondence 



between the highest weight states of the Yangian supermultiplets and the states of the CSM 
(i.e., they satisfy the same eigenvalue equations.) A further generalization to SU{n) case 
has also been accomplished in [|Ty,|rT|]. The lattice CSM model at even (odd) integer values 
of A can be mapped to bosonic (fermionic) spinless t-J model |[T^. 

29 



The method for evaluating the correlation functions of the CSM presented in this paper 
is directly applicable to its lattice cousins in some limited cases. The Galilean invariance is 
broken in the lattice models and is replaced with much weaker lattice translation symmetry 
which induces appearance of the Brillouin zone boundaries. If the elementary excitations 
created by a local operator acting on the ground state do not cross the Brillouin zone bound- 
aries and the excited states are the Yangian highest weight states, then the corresponding 
correlation functions of HSM are identical to that of the CSM. 



X. CONCLUSION 

In this paper the fractional exchange and exclusion statistics are studied using the exactly 
solvable Calogero-Sutherland model; and they are found to be mutually consistent. I have 
shown that the interaction giving rise to the fractional exclusion statistics for the elemen- 
tary excitations of a given condensed matter system can in fact be treated as the statistical 
gauge field carried by the real particles making up the system. This was done by calculating 
the exact dynamical density-density correlation function and one-particle Green's function 
using Jack symmetric polynomials and examining the intermediate states contributing to 
the correlation functions. I find that the intermediate states for X = p/q CSM are spanned 
by q quasiparticles and p quasiholes for the density-density correlation function and q — 1 
quasiparticles and p quasiholes for the hole propagator and thereby show that the quasipar- 
ticles indeed carry charge 1 and flux 2ttX and the quasiholes charge — 1/A and flux —2it as 
flrst suggested by Haldane 0. 

I also construct an explicit multi-sector anyon operators in analogy with Haldane's 
harmonic-fluid and calculate their correlation functions which agree with those of the 



CSM which corresponds to the flrst-quantized construction of anyons. Therefore, the CSM 
at odd-integer coupling constant describes the edge states of the fractional quantum Hall 
droplet corresponding to the Laughlin states as suggested in Ref. [ pT| . 
There are some interesting open problems: 
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• How can one rigorously construct ID fractional exchange statistics using the Yang- 
Baxter equation ? 

• For the CSM with the coupling constant other than A = 1/2 , 1, and 2, are there any 
corresponding random matrices? 

• How can one generalize Jack polynomials to the SU (n) case? 
I hope to see some of these questions answered in the future. 
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APPENDIX A: HOW TO TAKE THE THERMODYNAMIC LIMIT 

The thermodynamic limit of the DDDCF given by Eq. ( |5.3| ) greatly simplifies by the fact 
that the coefficient [0']^ vanishes unless the diagram 2^(k;) has at most p rows of length longer 
than q and q columns of length longer than p. As shown in Section jTV] the minimal charge 
neutral excitations in the CSM consist of q quasiparticles and p quasiholes; therefore, the 
intermediate states that contribute to the correlation function have exactly q quasiparticle 
and p quasihole excitations.^ 

The contributing diagrams T>{k,), for convenience, are divided into three subdiagrams 
A{k,) = l<t<p, l<j<q}, B{n) = j), 1 < J < g, g + 1 < ^ < and 

C(k) = {(i, j), 1<«<P, P+l<j< K^i} as illustrated in Fig. ^. The factors expressed 
in terms of the generalized factorial over 'D{n) are evaluated for each of the subdiagrams, 
separately. Then, later the each subfactors are multiplied to obtain to full factors over the 
diagram T). 



^ The non- vanishing diagrams with less than p rows and q columns are interpreted as having some 
of the quasiparticles and quasiholes in their unexcited modes. 
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First, consider the following factor that appears in the DDDCF 

J,Lv N+ji\-i ) 

Fi over A simplifies in the thermodynamic limit to 

In order to evaluate Fi over B I rewrite the product in terms of the gamma functions using 
the identity T{z + n)/T{z) = z{z + 1) ■ ■ ■ {z + n — 1) as 

^ A (iV+(j-l)/A-^;. + l)---(iV+(j-l)/A-g) 

_ r(iV + j/X - k'^)T{N + j/A - g + 1 - 1/A) 

Mr(iv + jyA-«:;. + i-i/A)r(iv + j7A-g)- ^ ^ 

Using the following asymptotic relation, 

hm ^ ^ = z"", (A4) 

IzHoo T[zj 

I reduce in Eq. as 

- n(l-«:;/iV)^/^-^ (A5) 

Fi over C, similarly, is given by 

9 / 1 X A-1 

i=i ^ 

Evaluation of ([0']^)^ is straightforward and is equal to product of the following three 
factors 

mknA) = A-(-i)r^(p) n plrl'llvA) ' ^^^^ 

- l\ P(-(,-l)/A + g)' 

= A-EL.'^O.... ft (A9) 



Fi(C)^^-nfl + i|) ~\ (A6) 
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The value of over A is given by 



JK{A) = f[f[{K'^ + \K,r. (AlO) 



J ■ \ ■ 



In order to evaluate over B, I further divide up the subdiagram B into p cells such 
that Ith cell is given by I < j < I, n'lj^i + 1 < ? < n'l}- If = /t^+i, then the /th cell 

is empty. Fig. ^ illustrates how the subdiagram B is divided into p cells and how the empty 
cells appear. The value of over B is then as follow 

j^(^) = nri n k-^+i+(^-j)/a)k-^+(^-j+i)a) 

l = lj = li = K[^^+l 

_ n n - + 1 + (/ - j)/m^', - + (/ - J + i)/A) 
y^^y r^-«:i + i + (/-j)/A)rK.-K; + (/-j + i)/A) ' ^ > 

where Kp^i = q- The contributions from the empty cells are identically equal to one. For 
the non-empty Zth cell Ki = I for + 1 < i < kJ. The expression above simplifies further 
to 

j^(^) = n n^j + 1 - j/m^- - u - i)/a) 

TjK', - + 1 - (z - J + l)/X)TiK[ - k;. - (z - j)/A) 

y r(/.^ + j)/A)r(/.'. - - - j - i)/a) ' ^ ^ 

In order to take the N ^ oo limit using Eq. ( [A4|) , I take the following ratio first and then 
send the ratio to the limit, 

Similarly, over C is given by 

= n ri'^. - A(^ - i))r(K, + 1 - AO 

^ y-^) i=i 

^ TT r(/t^ - /tj - A(i - j))T{Ki -Kj + 1- \{i - j + 1)) 

I take the following ratio and sent it to the thermodynamic limit using Eq. ([A4| ) , 

'^K^ic) N^oo r«(A) p 



U'^j-'Ul^^-'^A""- (A15) 

j=l i>j 
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Putting all the terms together, setting h^j/N = Xj and Kj/N = yj and turning the sums into 
integrals I get Eq. ( |5.4| ) with the normalization constant given by Eq. ( p.8| ). 

If the differences — — are of order 0{1) then the corresponding contribu- 

tions to Eq. ( |M3| ) (Eq. ( [XT5|) ) are of order C(l) instead of 0{N^/^) {0{N^^)) as iV ^ oo; 



therefore, the contributions of such terms to the DDDCF are suppressed in the thermody- 
namic limit. If n'j is of order 0{1) or 0{N) then there is a corresponding contribution of 
order 0{N^~^^^); so, if A < 1 the form factor vanishes as ^ oo, and if A > 1 it diverges. 
However, the net contribution of such term to the DDDCF is of order A^~^ x N^^^/'^ which 
always vanishes so long as A > 0^. (The factor comes from Ai/j.) More generally, let 
there be r k'j of all order 0{1) or 0{N). Then, the form factor will be of order N to the 
power of r(l — 1/A) — (2/A)r(r — l)/2 = r(l — r/A). Therefore, if the A > r the form factor 
diverges. In other words such configuration is favorable. This is an exotic violation of the 
super selection rule since r quasiholes are stuck together. However, the net contributions 
of these terms are again vanishingly small except, of course, in the long- wavelength limit 
where only these exceptional states survive. 

Almost identical method can be applied to obtain Eq. ( |6.5| ). When the limit x,t ^ is 
taken, the propagator becomes just the static density po- This gives a generalization of the 
Selberg's integral formula and might have connections with the g-deformed Lie algebra. 

APPENDIX B: ASYMPTOTIC EXPANSION 
1. Dynamical density-density correlation function 

A method for expanding the dynamical density- density correlation functions as x — > cxo 



is presented in the appendix. A similar method has been used previously in |]18[ for the 
equal-time correlation functions. 

The form factor to the DDDCF is largest near a phase space region where all Xj's near 
zero and y^'s near zero or one. Thus, the integrand will be expanded about this region. In 
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Fig. the shaded region gives non-zero contributions to the DDDCF and the regions with 
darker shades give largest contributions to the DDDCF in the long- wavelength limit. Each 
of these dark regions is labeled with the index n and is called nth sector, whose contribution 
to the correlation function has a characteristic oscillation, and is spanned by all Xj's close 
to zero and n {p — n) i/j's close to one (zero). 

To leading order in the normalized momentum variables the following relation hold true 



p—n 



Qx±Et = 2npoC \T.^i + T.yj] - 27rpoe^ E + 2vrpona;, (Bl) 
\i=i j=i J k=i 

where = x^Vgt. The sound velocity Vg is 27rpoA. A new set of n variables Wk for yj close 

to one are introduced such that Wk = 1 — yj- {^~) is the space-time coordinate for the 

right (left) moving particle or hole excitations. 

To leading order in each of the p sectors the DDDCF is given by 

(p(x,t)p(0,0)) ^ ^ (^I^Mq,p) +I^H\q,p) 

+ ECn (^^) ' I {l^i~t\q,p-n)I^i+t\n)e''''''^-^ 

+ Ii{^\q,P- n)/2^(-2|n)e-^2"^«'^^) , (B2) 

where C„ = (j^CX''^"''^n'^ (2^) ^ functions and I2 are defined as follow 



i=i 1 if m ^ 

2 Oi<j ^j\ Y[i'<j' 



X TT IT (3^^ + Xyj)-' ""'^ ^' /^-'<^''y^/^'' 

11 11 V « I yjj I 1-1/A 



xexp\-z\Y,x, + J2yA] (B3) 



I^{z\m) = n / dw, ^^'<!^ ' expi-zY^wA (B4) 

Since and I2 are absolutely convergent only if Re{z) > 0, it is necessary to analytically 
continue the functions as follow 
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/ 1 \ -^'^+'"^/-^~2im+25m p 

I^{z\l,m) = i^-j (B5) 
I^{z\m) = i^-j I^{l\m). (B6) 

The above analytical extensions give the following relations: Ii{—i\q,P — n)/2 = 
Ii{l\q,p- n)/2 (l|n) and Ii{±i\q,p) = -Ii{l\q,p). 

Finally, the DDDCF to leading order in each harmonic mode is given by 

c ( \ y''^^ 

(0|p(a;,t)p(0,0)|0) ^ f ^ + + X:c„ cos(27rpona;), 

(B7) 

where C„ = C„/^(l|g,p — 'n,)/2 (l|n). 

The renormalization parameter is defined in as r/ = 2(f j/wat)^/^, where wj and wat are 
the current and charge velocities. For the f/(l) CSM r] = 2/\ |lO|-0 With this identification 
the form of Eq. (|B7|) in the static limit agrees with the expansion given in Furthermore, 



by matching the universal constant in the 0th sector I deduce the following integral formula 

li{l\q,p) = {27i^\C)-\ (B8) 

2. One-particle Green's function 

Using the similar method used for the DDDCF I find the following leading order terms 
for the hole propagator 

(0|vl/t(x,t)vl/(0,0)|0)^5:i)„, e*(2-Po(n-A/2)x+(.poA)^t)^ ^gg) 

„=o \2vrpo4 + / / 

where Vn = poDX^-^-^^'^-^^'' (^l^I^{l\q - l,p - n)I^{l\n)e-'''^''-^^'^\ Here, I^{z\q - l,p) is 
defined as before but without the pre-factor + J2j Uj 



The interaction coupling constant A used in Ref. |1C| corresponds to A — 1 in this paper 
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FIGURES 

FIG. 1. A level diagram for a family with the highest-level mother-state given by 

|6)i = $(4,3, 1,0). The arrows comicct the mother states to their daughter states that are gener- 
ated by squeezing on the pairs of quantum numbers by one unit. 

FIG. 2. The shaded area is the region of support for the dynamical density-density correlation 
function at A = p/g = 5/3 in the energy- momentum space. The energy E is in some arbitrary unit 
and the momentum Q in units of 2ttpo. The 2p distinct regions with darker shade corresponds to 
the low-energy sectors of the model. 

FIG. 3. Cylinderical geometry for the anyon system. Two oppositely-moving edges along x-axis 
are separated by an incompressible anyon fluid. The cylinder is pierced by a solenoid which induces 
a magnectic flux $ through its cross-section. 

FIG. 4. A diagram V that contributes to the DDDCF at A = p/g = 5/6 is divided into three 
sub diagrams A, B, and C for convienience. 

FIG. 5. A subdiagram B with jo = 10 is divided into p cells. The dots • indicate the empty 
cells where k'i = and only the 1, 3, 4, 5, 8, and 10th cells are non-empty. For the Zth cell, Ki 
are all same and equal to I. 
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